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ABSTRACT 

Context. M dwarfs have been found to often have super-Earth planets with short orbital periods. Such stars are thus preferential targets 
in searches for rocky or ocean planets in the solar neighbourhood. 

■ Aims. In a recent paper (Bonfils et al. 201 1), we announced the discovery of respectively 1 and 2 low mass planets around the M1.5V 
[•-p'] I stars Gl 433 and Gl 667C. We found those planets with the HARPS spectrograph on the ESO 3.6-m telescope at La Silla Observatory, 

, from observations obtained during the Guaranteed Time Observing program of that instrument. 

• Methods. We have obtained additional HARPS observations of those two stars, for a total of respectively 67 and 179 Radial Velocity 
CLf measurements for Gl 433 and Gl 667C, and present here an orbital analysis of those extended data sets and our main conclusion about 

both planetary systems. 

Results. One of the three planets, Gl 667Cc, has a mass of only M2. sin ; ~ 4.25 and orbits in the central habitable zone of its host 
star. It receives just 10% less stellar energy from Gl 667C than the Earth receives from the Sun. However planet evolution in habitable 
^ ■ zone can be very different if the host star is a M dwarf or a solar-like star, without necessarily questioning the presence of water. The 

two other planets, Gl 433b and Gl 667Cb, both have M2. sin i of ~5.5 Me and periods of ~7 days. The Radial Velocity measurements 
of both stars contain longer time scale signals, which we fit as longer period Keplerians. For Gl 433 that signal probably originates 
in a Magnetic Cycle, while a longer time span will be needed to conclude for Gl 667C. The metallicity of Gl 433 is close to solar, 
^ ■ while Gl 667C is metal poor with [Fe/H] — 0.6. This reinforces the recent conclusion that the occurence of Super-Earth planets does 

' not strongly correlate with stellar metallicity. 

■ Key words, techniques : radial velocity / stars : late-type / planetary systems 
'NT ■ 

(N 

. 1 . Introduction P<10 days, to be compared to ~ 50% for similar planets around 



o 



■ Much interest has recently focused on planets around M dwarfs. 



G dwarfs. 



X 



. , , . . . . . ^ . , . The physical conditions in the circumstellar disks of very 

with three mam motivations: constraining formation physico- ^^^^ ^^^^ jj^^^^f^^^ ^^^^^ ^j^^ formation of low-mass planets 

chemical characterization of planets and finding rocky planets ^^^^j^ ^^^^ ^^^^^ ^^^^^^^ ^1^^^ ^j^^ ^^^^ ^j^^^ translates 

m the habitable zone of their stars The compared occurence fre- .^^^ ^^^^ ^^jj^^^^ ^^^^^^ transiting M dwarfs 



quency, as a function of orbital elements, of planets around M ^j-j^^ ^j^^^^ ^.^^^^^^ announced by Muirhead et all 12011 . 

; dwarfs and around the more massive solar-type stars probes the ^^^j^ transiting planets ai-e excellent atmospheric characteriza- 

■ ■ sensitivity of planetary formation to its i nitial conditio ns^t ^.^^ ^ ^.^^^ ^^.^^^.^ ^ ^^^j^^ ^.^^ .^^ ^^^jj^^. 

planets are rare around M dwarfs, with j Bonfils et al. | miM), radius, a transit across an M dwarf provides much more accurate 

"'^^'J^"'^'"^ \ frequency of 6! % for periods under i^easurements of radius and transmission spectrum than across 

10000 days. That number IS lower than the 10+2% frequency for , 

1 111-1 / rrr — — a ^r\i 1 a soiar-iype star, 
similar planets around solar-like stars (e.g. IMavor et al. 2011), 

though not yet at a high significance level. Super-Earth (2- Finding rocky planets in the habitable zone (HZ) of their 

lOMe), by contrast, seem abundant around M dwarfs at the short ^tars is another motivation for planet searches around M dwarfs, 

orbi tal periods to which ra dial velocity searches ai-e most sensi- Pl^^ts of given mass and orbital separation induce larger stel- 

tive: lBonfilsetal.1 SMU find an occurence rate of 88+^^% ^r radial velocity variation around lower mass stars, but, more 

importantly, the low luminosity of M dwarfs moves their habit- 
able zone much nearer to the star These two effects combine. 



* Based on observations collected with HARPS instrument on and a habitable planet around a O.S-Mq M dwaif produces a 
the 3.6-m telescope at La Silla Observatory (European Southern 7 times larger radial velocity wobble than the same planet or- 
Observatory) under programs ID072.C-0488(E) biting a solar-mass G dwarf. Additionally, some atmospheric 
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models suggest additional advantages o f planets in the HZ of M 
dwarfs for habitability characterization: Segura et al. ("2005") find 
that the significant dependence of atmospheric photochemistry 
on the incoming spectral energy distribution strongly reinforces 
some biomarkers in the spectra of Earth-twin planets orbiting M 
dwarfs. N2O and CH3CI, in particular, would be detectable for 
Earth-twin planets around M dwarfs but not around solar-type 

stars. 

As discussed in iBonfils et al.l (1201 lal) . which presents the 
full specifications of the survey, we have been monitoring the 
radial velocities (RV) of a distance and magnitude limited sam- 
ple of 102 M dwarfs since 2003 with the HARPS spectrograph 
mounted on the ESO/3.6-m La Silla telescope. With a typical RV 
accuracy of l-3m/s and 460 hours of observations, our survey 
identified super-Earth and Neptune-like planets around Gl 176 
(Forveille et al. 2009), Gl 581 (Bonfils et al. 2005b; U drv et alJ 
2007; Mavoretal. 2009) and Gl 674 (B onfils et al.li2007h . In 
Bonfils et al., (.20 11 a) , which was centered on the statistical im- 
plications of the survey for planet populations, we added to this 
list 1 planet around Gl 433 and 2 around Gl 667C, providing de- 
tailed periodograms and a false alarm pro bability in using our 
GTO/HARPS data. lAn^lada-Escude et all (1201 2i) recendy an- 
nounced a confirmation of the planets Gl 667Cb and Gl 667Cc, 
which howe ver is only very partly independent, since it largely 
rests on our lBonfils et al.l (1201 lah observations and data reduc- 
tion. 

Here we present a more detail ed orbital analysis of t hese two 
systems than we could present in lBonfils et al.l (1201 lah . and we 
refine their parameters by adding new seasons of RV measure- 
ments. The 3 super-Earths have M2. sin / between 4.25 and 5.8 
Me. With periods of ~7days, both G1433b and G1667Cb are hot 
super-Earths. G1667Cc, by contrast, has a 28-day period. It or- 
bits in the centre of the habitable zone of its star, and receives just 
10% less stellar energy than the Earth receives from the Sun. We 
also detect some longer-period radial-velocity variations, which 
we also discuss. 

The next section discusses our data taking and analysis, 
while Sect. 3 summarizes the stellar characteristics of Gl 433 
and Gl 667C. Sect. 4 and 5 respectively present our orbital anal- 
yses of the planetary systems of Gl 433 and Gl 667C. Sect. 6 
discusses the habitable zones of M dwarfs, with emphasis on the 
case of Gl 667Cc. Finally, Sect.7. summarizes our conclusions. 

2. Spectroscopic and Doppler measurement with 
HARPS 

Our observin g proc edure is presented in some details in 
IBonfils et all (l2011ah . and is only summarized here. For both 
stars we obtained 15 min exposures with the HARPS spec- 
trograph (High A ccuracy Radial velocity Planets Searcher 
iMavoretaT] 12003b . HARPS is a fixed-format echelle spectro- 
graph, which covers the 380 to 630 nm spectral range with a re- 
solving power of 1 15 000. HARPS is fed by a pair of fibres, and 
is optimized for high accuracy radial-velocity measurements, 
with a stability better than Im/s during one night. To avoid light 
pollution on the stellar spectrum of our "faint" M-dwarf targets, 
we chose to keep dark the calibration fiber of the spectrograph. 
Our Radial Velocity accuracy is therefore intrinsically limited 
by the instrumental stabiliy of HARPS. That stability is however 
excellent and the signal-to-noise ratio of our 15min exposures 
of these two V~10 M dwarfs limits the RV precision to slightly 
above 1 m/s (the median S/N ratio per pixel at 550nm is 57 and 
65 respectively for Gl 433 and Gl 667C). Sect. 4 and 5 discuss 
the number of exposures and their time span for the two stars. 




Fig. 1. Emission reversal in the Ca II H line in the average 
spectra of (from the less to the more active) Gl 581 (M3, in 
magenta), Gl 667C (Ml. 5, in pink), Gl 433 (Ml. 5, in yellow), 
Gl 176(M2.5, dark) and Gl 674 (M3, in blue). Within our 100 
M dwarfs sample, Gl 581 has one of the weakest Ca II emission 
and illustrates a very quiet M dwarf. Gl 674 and Gl 176 have 
a much stronger emission and are both moderately active with 
an identified rotational period of respec tively 35 and 39 days 
(IBonfils et al.ll2007l:lForveille et al.ll2009l) . 



For homogeneity, we reprocessed all spectra with the lat- 
est version of th e standard HARPS pipeline. The pipeline 
(iLovis & Pepe 2007) uses a nightly set of calibration exposures 
to locate the orders, flat-field the spectra (Tungsten lamp illu- 
mination) and to precisely determine the wavelength calibration 
scale (ThAr lamp exposure). We measured the Radial Velocity 
through cross correlation of the stellar spectra_with a numerical 
weighted mask, following the procedure of jPepe et al.l (l2002h . 
For both Gl 433 and G1667C, we used a mask derived from a 
very high S/N spectrum of a M2 dwarf. The velocities of Gl 433 
and Gl 667C have internal median errors of respectively 1.15 and 
1.30 m/s. This includes the uncertainty of the nightly zero point 
calibration, the drift and jitter of the wavelength scale during a 
night, and the photon noise (the dominant term here). 

Because both stars have significant proper motion, the pro- 
jection of their velocity vector changes over the d uration of our 
survey. We subtract this secular acceleration (see iKiirster et al.l 
2003, for details) before radial velocity analysis. For Gl 433 and 
Gl 667C the values are respectively 0.15 and 0.21 m/s/yr. 



3. Stellar characteristic of Gl 433 and Gl 667C 

Both Gl 433 (LHS 2429) and Gl 667C are early-M dwarfs in 
the close solar neighborhood. Table [T] summarizes their proper- 
ties. The masses are computed from the /T-band absolute mag- 
nitudes using the jPelfosse et al.l (2000) empirical near-infrared 
mass-luminosity relation. The bolometric correction of Gl 433, 
and then its luminosity, is computed from the / - K color us- 
ing the cubic polynomial of iLeggett et al.l (12000). Those authors 
directly determined the luminosity of Gl 667C using a combi- 
nation of flux calibrated observed spectra and synthetic spectra, 
and we adopt their value. 
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Spectral type 
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T(2) 
J 




f. s5+n 09 

U.OJItU.UZ. 
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J.OUIIIU.UH 




1^(2) 
JV 


J.DZltU.UZ 


A n^-t-o 07 

D.UJItU.UZ 


71 


110.6±1.8'^^' 


138.2+0.7*'" 


M (in Mo)*'' 


0.48 


0.33 


L (in Lo)^' 


0.034 


0.014 


[Fe/H]<''> 


-0.22 


-0.55 


TefF'*) 


3600 


3600 


logi^x'" 


<-4.8 


<-4.12 



Table 1. (1) Hawley et al. ( 199 1); (2) ICutri etal.1 (l2003h : (3) 
iPerrvman et all (Il997i) : (4) Soderhjelm (1999*); (5) from mass- 



lumino sity relation of iDelfo sse et al. (2000) ; (6) Leggett et al 
20001) : (7) from luminosity-c olor relat ion of Bonfils et al 



J2005ah : (8)|Morales elaD (l200 8j); (9) fromlSchmitt et alj (i 19951) 
for Gl 433 and from Schmitt & LieflS (l2004|) for Gl 667C 



3. 7. G/ 433 ; metallicity, activity and dynamic population 

According to iBonfils et al.l (l2005ah we estimated the metallic- 
ity of Gl 433 to [Fe/H] — 0.2. This photometric calibration rela- 
tionship is based on a faint number of stars in the high metal- 
licity range but Neves et al. (2012) show that for such value 
of [Fe/H] the calibration is correct with a typical dispersion of 
0.2 dex. This is confirmed by the value of [Fe/H] =-0. 13 obtained 
in us ing Neves et al. (2012) relationship, itself an update from 
ISchlaufman & LaughlinI (12010). We could conclude that Gl 433 
is solar or slightly sub-metallic. 

Gl 433 is d ynamically classified as a membership of old disk 
(lLeggettlll992|). Ca II H and K chromospheric emission is de- 
termined bylR^scheT&Marc^C 2006) and Gl 433 is in the less 
active half of M dwarfs with same luminosity. For objects of 
similar spectral type a direct comparison of Ca II emission lines 
gives relative estimation of rotational period (larger Ca II emis- 
sions correspond to shorter rotational periods). In the Fig [1] we 
present an average HARPS spectrum of Gl 433 in the region of 
Ca II H line. Gl 433 has an emission of Ca II slightly higher 
than Gl 581 (a very quite M dwarf), but well below those of 
Gl 176 and Gl 674 for whose rotational pe riods of respectively 
35 and 39 d ays have bee n determined by iBonfils et all (120071) 
and Foi-veille etaLl (l2009h . The five M dwarfs of the Fig. [1] do 
not have identical spectral (from Ml. 5 to M3) therefore estima- 
tion of rotation from the Ca II emission is only indicative. It 
shows that Gl 433 rotates most probably with a period longer 
than Gil 76 and Gl 674. 

The X-ray flux of Gl 433 is not detected by ROSAT and 
we use a ROSAT limiting sensitivity of 2.5^^ .[d/ 10 pc\^ erg/s 
(ISchmitt et al.l[T995i) to estimate Rx = LogLxlLsoL < -4.8. For 
an M dwarfs of ~ 0.5Mc7) th e Rx versus rotation period relation 
of iKiraga & Stepieiil (l2007l) give Proi > 40 days for such level 
of X flux. This inferior limit on the rotation period is coherent 
with the estimate from the Ca II emission. 



3.2. Gl 667C : multiplicity, metallicity, activity and dynamic 
population 

Gl 667C is the lightest and isolated component of a hierarchical 
triple system, the two others components are a closest couple of 
K dwarfs. Gl 667 AB has a semi-major axis of 1.82 A.U. (pe- 
riod of 42.1 5 years) and a total mass dynamically determined 
of 1.27Mo (ISoderhielmlll999i) . Gl 667C is at a projected dis- 



tance of 32.4" of Gl 667 AB, giving an expected semi-major axis 
of ~300 A.U. (for a distance of 7.23 pc and a factor of 1.26 be- 
tween expected and projected semi-major axis, Fischer & Marc^ 
Il992h . 

The IBonfils et al.l (12005 ah photometric relationship for 
G1667C gives us an estimate of the metallicity of 0.55 
dex, which agrees quite well with the spectroscopic deter- 
mination of [Fe/H]~-0.6 for the primary (Perrin et al. 198.^ 
Zakhozhai & Shaparenko 1996). This demonstrates that the 
Bonfils et al., (,2005 a) photometric calib ration gives excelle nt re- 



sults for low-metallicity M dwarfs. The lNeves et al.l (1201 2h rela- 
tionship gives [Fe/H] — 0.45, confirming that this star is a metal- 
poor M dwarfs. 

Consistent with this chemical composition, G1667 is classi- 
fied as member of the old disk population from its UVW veloc- 
ity iLeggett 1992) and is among the objects of lower chromo - 
spheric emission for its luminosity (Raus cher & Marcvll2006l) . 
In the Fig. [U the Ca II H line emission of Gl 667C is slightly 
inferior to the Gl 433 one, also indicating also a large rotational 
period . 

The NEXXUS database (ISchmitt & Liefkdllool indicates 
detection of the X-ray emission of Gl 667C from the ROSAT 
All-Sky Survey (RASS) with a value of log Lx = 27.89. But the 
resolution of ROSAT image is not sufficient to separate the three 
components of the system. Although this value is the jointed flux 
for Gl 667 ABC, it gives a superior limit on the coronal emis- 
sion of Gl 667C. This implies Rx - LogLx/LsoL well below 
-4.12 and rota tional period superior to 40 days (in applying the 
iKiraga & Stepi en (2007) Rx versus rotation period relationship). 

4. Orbital analysis of Gl 433 

4.1. HARPS measurements 

We obtained 67 measurements of G1433 radial velocity that span 
2904 days from December 2003 to November 201 1. The data 
add 17 new points to those analyzed in lBonfils et al.l(l201 lab and 
extend their time span by 1200 days. We estimate that the uncer- 
tainties of ~1.1 m/s are dominated by photon noise (the median 
S/N ratio per pixel at 550 nm is 57). Overall, the RVs have a 
rms=3.14 m/s and - 7.7 + 0.5 per degree of freedom, con- 
firming variability in excess to measurement uncertainties (~1.1 
m/s, dominated by photon noise). 

We therefore continued with a classic al periodicity analy- 
sis, based on floating-mean p eriodograms (Gilhland & BaliunasI 
ll987l:IZechmeister et al.ll2009i) . Fig.|2](bottom panel) depicts the 
periodogram of our RV time series and shows a strong power 
excess around the period 7.3 day. We plot the window function 
in Fig. |2] (middle panel) to identify the typical time sampling but 
found no counter part to this periodici ty. We adopted the normal - 
ization of periodograms proposed by Zechmeister et alj ( l2009l) . 
which is such that a power of 1 means that a sine function is a 
perfect fit to the data, whereas a power of indicates no improve- 
ment over a constant model. Hence, the most powerful peak was 
measured at a period P~7.3 days with a power p,„a.x - 0.51. 
To assess its fortuity we performed a bootstrap randomization 
(Press etal. 1992): we generated 10,000 virtual data sets by 
shuffling the actual radial velocities and retaining the dates: for 
each set we computed a periodogram: with all periodograms we 
built a distribution of power maxima. It appears that, in ran- 
dom data sets, power maxima exceed p = 0.40 only once every 
100 trials, and never exceed p - 0.43 over 10,000 trials. This 
suggests a False Alarm Probability (FAP) < 1/10,000 for the 
7. 3-d period seen in the original periodogram. Moreover, for a 
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Fig. 2. Top : HARPS radial velocity of Gl 433. Middle : win- 
dow function of the measurement. Bottom : Periodogram of the 
HARPS measurements. 

more prec i se est imate of this low FAP value, we made use of 
ICummingl(l2004 analytic formula : FAP = M.(l - p„,axf^'^^^^, 
where M is the number of independent frequencies in the peri- 
odogram, p,„ax its highest power value and the number of mea- 
surements. We approximated M by 2904/1 (the ratio between the 
time span of our observations and the typical 1-day sampling), 
and obtained the very low FAP value of 4 x 10"^. 

We modeled that first periodic signal by a Keplerian orbit. 
Even though the periodicity was securely identified we used 
Yorbit (an heuristic algorithm, mixing standard non-linear min- 
imisations and ge netic algorithms; Segransan et al. in prep; 
iBonfils et al.ll201 fa ) and benefited from a global search without 
a priori. We converged on a period P - 7.3732 ±0.0023 d, semi- 
ampUtude Ky - 2.99 + 0.38 m/s and eccentricity e = 0.17 + 0.13 
which, for M^, - 0.48 Mq converts to m sin / = 5.49 + 0.70 Me 



(Table |2|. The rms and ^Jx^ around the solution decreased to 
2. 17 m/s and 2.00 + 0.09 per degree of freedom, respectively. 

Inspecting the residuals and their periodogram, we measured 
the power of the most powerful peak {P ~2900 d; - 0.28) 
and apply bootstrap randomization to find it is insignificant 
{FAP - 22%). We nevertheless found that adding a quadratic 
drift to the 1 -planet model did improve the solution (rms=1.91 
m/s and ^[x^ - 1 .79 + 0.09 per degree of freedom). 

4.2. HARPS+UVES measurements 

Radial vel ocity measurement s of Gl 433 were also obtained with 
UVES by IZechmeister et al.1 (l2009l) . They span 2553 days be- 
tween March 2000 and March 2007 and can therefore extend the 
time span of our observations by more than 3 years (see Fig.[3]- 
second panel). 

Hence we pooled together UVES and HARPS data and 
performed again a 7 planet fit with Yorbit. We converged on 
the same solution with more precise orbital parameters (P - 
7.37131 + 0.00096 d; Ki = 2.846 + 0.248 m/s; e = 0.13 + 0.09). 
The residuals around the solution show a power excess at a sim- 
ilar period than before {P ~2857 d; p,„ax - 0.19) that now ap- 
pears more significant. With 10,000 trials of bootstrap random- 
ization we found a FAP ~ 0.2%. Although such signal requires 
additional measurements for confirmation it justified the addi- 
tion of a second planet to our model. 

With a model composed of 2 planets on Keplerian orbits, 
Yorbit converges on periods of 7.37029+0.00084 and 3690+250 
days, with semi-amphtudes of 3.113+0.224 and 3.056+0.433 
m/s and eccentricities 0.08277+0.07509 and 0.17007+0.09428 
(see. Table [3] for the full set of parameters). The solution has a 
rms=2.39 m/s and a ^[x^ - 1.30 + 0.05 per degree of freedom, 
showing that the combination of HARPS and UVES data sets 
largely improved the parameter uncertainties. 

4.3. The planetary system around Gl 433 

In the past, Gl 433 was announced as hosting a brown dwarf of 
30 Jupiter mass with an orbital p eriod of ~500 days from astro- 
metric measurements iBernsteinI (Il997h . But all radial velocity 
measurments published after, as the present one, reject this de- 
tection. 

Apparent Doppler shifts may also originate from stellar sur- 
face inhomogeneities, such as plages and spots, which can break 
the balance between light emitted in the red-shifted and the 
blue-shifted parts of a rotating star (e.g. lSaar & Donahud[T997l: 
[Oueloz et al. 2001; Desortetal. 20075. However the large ro- 
tation period (>40 days) of Gl 433 ensures that the observed 
doppler shift of 7.37 days does not originate from the stellar 
activity. A search of correlation between the 7.37 days radial 
velocity period and H„, Call-index or bissector is unsuccess- 
ful in our HARPS data. Such signal is easily detected when 
the radial velocity variation is due to activity fo r a period of 
~35 days (an d a fortiori f or all shorter periods) teonfils et al.l 
2007; Forveil le et al]|2009l see the case of Gl 176 and Gl 674). 
This is a strong evidence that the ~7 d period is not due to stellar 
surface inhomogeneities but to the presence of a planet orbiting 
around G1433. 

This planet, Gl 433b, belongs to the category of super-Earth 
with a minimum mass of ~5.8 Mg. At a separation of 0.058AU 
from its star, Gl 433b is illuminated by a bolometric flux, per 
surface unity, 10 times higher than what the Earth receives. 
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Table 2. Fitted orbital solution for the HARPS data of Gl 433 
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Fig. 3. 2-planets model of HA RPS (in red) and UVES (in blue, 
from lZechmeister et al.ll20()P ) radial velocity measurements for 
Gl 433. Top : phased radial velocity for a planet of 7.37 days 
period. Middle : radial velocity curve for the 3700 day period. 
Below : residuals and their periodogram. 
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Fig. 4. Top : HARPS radial velocity of Gl 667C. Bottom : one 
planet plus a linear drift model. 



The long period variability detected in the HARPS and 
UVES data set may have several origins. The signal can have 
as origin a ~10-year period planet of 50 M®. A possible ori- 
gin may also be an effect of a long term stellar magnetic cycle 
(like the so-called 11 -year solar cycle) during which the frac- 
tion of the stellar disk covered by plages varies. In such re- 
gion the magnetic field attenuates the convective flux of the 
blueshifted plasma and impacts the mean observed radial ve- 
locity in the order of few m /s or more (iMeunier et al.li2010l) . 
iGomes da Silva et al.l (|201 ih do detect large period variation for 
numerous activity proxy (Call, Nal, Hel, Hq, lines) of Gl 433 
coherent with a period of ~10 years. The RV variability at long 
period is well correlated with these activity proxies (Gomes da 
Silva et al. in prep). Therefore, we favor a Magnetic Cycle as the 
origin of the long term signal. 
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Table 3. Fitted orbital solution for the HARPS and UVES data Table 4. 1 planet + 1 linear drift orbital solution for Gl 667 
ofGl 433 
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M2. sin; [Mg,] 


5.79 


44.6 
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0.0504 


a[AU] 
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3.6 


y [km/s] 


6.55+0.02 


N 
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7 [m/s/yr] 


1.26+0.18 


S pan [days] 
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r.m.s [m/s] 


2.51 


r.m.s. [m/s] 

1 T 
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V? 
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5. Orbital analysis of Gl 667C 

We obtained 179 measurements of Gl 667C radial velocity span- 
ning 2657 days between June 2004 and September 2011. The 
data adds 36 points to those analyzed in iBonfils et alJ (1201 lah 
and extend their time span by 1070 days. They have uncertainties 
of ~ 1 .3 m/s (dominated by photon noise) and both rms and ^Jx^ 
values (resp. 4.3 m/s and 3.04 per degree of freedom) indicate 
a variability above those uncertainties. All RVs are shown as a 
function of time in Fig.|4](top panel) whereas their window func- 
tion and floating-mean periodogram are presented Fig. |5] (top 
two panels). 

5.1. A one planet and linear fit solution 

The long-term drift seen in Gl 667C RVs agrees well with the 
line-of-sight acceleration induced by its companion stellar pair 
Gl 667 AB which is about GMabI r\g_(^ ~ 3m/s/yr (for a to- 
tal mass Mab of 1.27Mo and a separation of ~300 AU). Even 
without removing that linear drift the periodogram additionally 
shows a peak at P~7.2 days (see Fig. |5] panel b). We do re- 
move an adjusted drift (7 = 1.59 + 0.14 m/s/yr) to make the 
periodic signal even stronger (pmax = 0.53; Fig.|5]panel c). To 
measure the FAP of the 7.2 days signal we ran 10,000 bootstrap 
randomization and found no power stronger than 0.18 suggesting 
aFAP « 1/10, 000. We also computed the FAP with jCummingl 
l2004 's prescription: FAP = M.(l - p.nax)*'^"^'^^- We approxi- 
mate M by 2657/1 (the ratio between the time span and the typ- 
ical sampling of our observations) and obtained the extremely 
low FAP value of ~ 10"^^. Note that, on Fig.|4](panel c), a sig- 
nificant and much less powerful peak is seen around P= 1 .0094 
day and corresponds to an alias of the 7.2-day peak with the typ- 
ical ~1 day sampling. 

We pursued by adjusting with Yorbit the RVs to a model 
composed of 1 planet plus a linear drift and converged robustly 
on the orbital elements reported in Table. |5] This model reduced 
the rms and ^jjfl per degree of freedom to 2.51 m/s and 2.02, 
respectively. They nevertheless remain above the photon noise 
and instrumental uncertainties, what prompted us to continue the 
analysis with the residuals and try more complex models. 

5.2. A multi-keplerian plus a linear fit solution 

Following up on the periodogram of the residuals around the 
7 planet + drift model (Fig |5] panel d), we observed additional 
power excess around a series of periods, the 6 most power- 
ful peaks being around 28, 91, 105, 122, 185 and 364 d, with 
p =0.225, 0.215, 0.185, 0.181, 0.134 and 0.123, respectively. 



Bootstrap randomization indicates FAPs lower than 1/10,000 for 
the 4 most powerful peaks whereas Cumming's prescription at- 
tributes them a FAP< 10"'^. 

To model the RVs with 2 planets + 1 drift, one could pick 
the most powerful peak as a guessed period and perform a lo- 
cal minimization to derive all orbital parameters. We found it is 
however not the most appropriate approach here because some 
of the less powerful peaks actually correspond to a signal with 
high eccentricities, and eventually turned to be better fit. To ex- 
plore the best solutions we generally prefer the global search 
implemented in Yorbit with eccentricities left to vary freely. We 
identified several solutions with similar ^[x2 values, reported in 
Table|5] For all solutions, the first planet and the linear drift keep 
similar parameter values (see Sect 5.1). On the other hand, the 
different solutions are discriminated by different orbital periods 
for the second planet, with P,. equals to -28, 90, 106, 124, 186 
or 372 days. 

To understand those signals we look at the periodograms of 
the residuals around each solution. We found that, on the one 
hand, when the solution is with Pc - 28 d, the periodogram of 
the residual shows the remaining peaks at 91, 105, 122, 185 and 
364 days (Fig.|6]- panel a). On the other hand, when the solution 
is with any of the other periods Pc = 90, 106, 124, 186 or 372 d 
then, only the 28 d peak remains (e.g. Fig.|6]- panels b, c, d, e and 
f). This means that all peaks at 90, 106, 120, 180 and 364 day 
actually correspond to a single signal, with its harmonics and 
aliases, and that the peak around 28 d is another independent 
signal. 

One of the 4 signals identified so far (3 periodic signals + 1 
linear drift) have ambiguous solutions. Among the possible pe- 
riods, the ~106-d period seems to correspond to the rotational 
period of the star as seen in one activity indicator (see. Sect. 
5.3). Conservatively, we continue by assuming that signal is due 
to activity and used the 106-d period to present our fiducial solu- 
tion. That assumption is given further credit when we run Yorbit. 
Without any a priori on any signal, Yorbit converges on a solu- 
tion with P - 1.2, 28 and 106 d. We present the orbital parame- 
ters of a 5 planet + 1 drift model in Tableland the RV decom- 
position and periodogram of the residuals around that solution 
in FiglT] That solution has rms and equal to 1.73 m/s and 
1 .44+0.06 per degree of freedom, respectively. 

Finally, we inspect the residuals around that last model (see 
panels (e) and (f) in Fig. |7]i. We found the power maximum of 
the periodogram located around P - 1.0083 d (a possible 1-day 
alias with a period of 121 day) with power /7,„f„ = 0.13, to which 
we attributed a FAP of 2.8%. Therefore, we did not consider that 
significant sine signal remains in the data. Nevertheless, remain- 
ing power excess around the period ~90 and 122 d suggests the 
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Fig. 5. From the top to the bottom : (a) Window function of the RV measurement; (b)periodograms of the Gl 667C RV measure- 
ments; (c) periodograms of the Gl 667C RV measurements after a drift removal; (d) periodograms of the residual after subtraction 
of a drift + 1 planet model. 



Table 6. 3 planets + 1 linear drift orbital solution 
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1.73 
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keplerian fit does not fully account for the signal we attributed 
to stellar activity. 

5.3. The planetary system around Gl 667C 

To assess if one of the radial velocity periodic signal could be 
due to stellar rotation we searched for periodicity on several 
activity diagnostics (Ha and Call-index, bisector-inverse slope 
(BIS) and full-width at half-maximum (FWHM) of the cross- 
correlation function). The clearest signal detected is a very high 
peak in the periodogram of the FWHM at a period of ~ 105 days 
with a FAP lower than 0. 1 % (see Fig. |8]l. 

P eriodic variation of the FWHIVI is an estimator of the rota- 
tion (Oueloz et al.ll2009l) . A rotational period of ~105 days for 
Gl 667C is consistent with the faint activity level of these stars, 
discussed in Sect. 3.2 and such period matches one of the power 
excess seen in our radial velocity periodogram analysis. Thus 



our favorite explanation is that the period of ~105 days and all 
the harmonics or alias at f = 90, 124, 186 and 372 days (see 
Sect. 5.2) are due to stellar rotation. However, we tried to use 
different sub-sets of RV points and found that the best solution 
may appear with Pj - 90, 106 or 186 d. We therefore consider 
that the RV signal at large (P > 90) period is not definitively as- 
signed to date, its origin may be stellar rotation (the most proba- 
ble in our point of view), or a high eccentricity planet, 2 planets 
in resonnance or a combination of several of these explanations. 

The RV signal at 7.2 and 28.1 d are completely independent 
of the large period ones and cannot be due to stellar rotation. 
Thus at least 2 super-Earths are present at close separation of 
Gl 667C wit h a mass of 5.5 and 4. 25 M^. Akeady announced by 
our team in lBonfils etaP (1201 id) with GTO-HARPS data, the 
use of supplementary HARPS RV allows us to specify their pa- 
rameters. At a separation of 0.05 and 0.12 a.u. from their star, 
Gl 667Cb and Gl 667Cc are respectively illuminated by a bolo- 
metric flux, per surface unity, 5.51 and 0.89 times than what 
the Earth receives from the Sun (in using stellar luminosity of 
Table [T]|. The super-Earth Gl 667Cc is then in the middle of the 
habitable zone of its star. We discuss this point in detail in the 
next section. 

6. A planet in the middle of the habitable zone of an 
M dwarf 

6.1. GI667CC 

Super-Earths in the habitable zone of their host stars, which 
by definition is the region where liquid water can be stable 
on th e surface of a rocky planet (Huang 1959; Kasti ng et al.l 
Il993h . currently garner considerable interest . For a detai l ed dis- 
cussion of the HZ we refer the reader to ISelsis et alJ (l2007h 
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Fig. 6. Periodograms for the Gl 667C residual of radial velocity measurements after subtraction of a drift plus a 2 planet solution. 
For all the panels the period of the first planet is Pt, - 7.2 days. From the left to the right and from the top to the bottom the period 
of the second planet is = 28, 90, 106, 124, 186 and 372 days. 



or lKaltenegger & Sasselovl (1201 l l). but we summarize the most 
salient points. One important consideration is that a planets with 
masses outside the 0.5-lOMe range cannot host liquid surface 
water. Planets under the lower end of this range have too weak 
a gravity to retain a sufficiently dense atmosphere, and those 
above the upper end accrete a massive He-H envelope. In ei- 
ther case, the pressure at the surface is incompatible with liquid 
water. The lOMg upper limit is somewhat fuzzy, since planets in 
the 3-lOMe range can have very different densities (reflecting 
different structure) for a given mass: Earth-like, Neptune-like, 
and Ocean planets can all exist for the same mass (e.g. Fig. 3 in 
[Winn etal. 2011). 

To potentially harbor liquid water, a planet with a 
dense atmosphere like the Earth need s an equi l ibrium 
temperature between 175K and 270K (ISelsis et al.1 l2007t 
[Kaltenegger & Sasselov 201 1, and references therein). If T^g > 
270K, a planet with a water-rich atmosp here will experien ce a 
phase of runaway greenhouse effect (see iSelsis et al.ll2007h . 



the other end of the range, CO2 will iiTeversibly freeze out from 
the atmosphere of planets with T^g < 175K, preventing a suf- 
ficient greenhouse effect to avoid freezing of all surface water. 
Either case obviously makes the planet uninhabitable. The inner 
limit of the habitable zone is well constrained. The outer one, by 
contrast, is very sensitive to the complex and poorly constrained 
meteorology of C02 clouds, through the balance between their 
reflecting efficiency (which coo ls the planet) and their green- 
house effect ( which warms it) (iForget & Pierrehumbertl [19971 : 



iMischna et al.ll2000l:ISelsis et al.ll20b7h . Finally, a location inside 
the habitable zone is a necessary condition for hosting surface 
water, but by no means a sufficient one: long term survival of 
surface water involves complex ingredients such as a carbonate- 
silicate geological cycle, and an adequate initial H2O supply. 

Adopting the notations of lKaltenepger et al] (1201 lb the equi- 
librium temperature of a planet is: 



At r,, = ((l-A)L„„,/(4ySD^)) 



.2^^1/4 



(1) 
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Fig. 7. Decomposition of our nominal 3 planet + 1 drift model for Gl 667C. (a) panel shows the contribution of the long-term drift 
with the 3 Keplerian contribution removed. Panels b, c and d show respectively radial velocity curves for keplerian f /, - 7.2 d, 
Pc = 28.1 d and - 106.4 d. Residual of this solution is showed in panel e with their periodogram in panel f. 



where ji is the geometrical fraction of the planet surface 
which re-radiates the absorbed flux (/S = 1 for a rapidly ro- 
tating planet with a dense atmosphere, like the Earth; fi - 0.5 
for an atmosphere-less planet that always presents the same 
side to its star) and A the wavelength-integrated Bond albedo. 
Using the stellar parameters of Table [1] the equilibrium temper- 
ature of G1433b and G1667Cb are respectively 495(^)'^'*K and 

426(^)'^'*K. Both planets are firmly in the hot super-Earth cat- 



egory. Gl 667Cc, by contrast, has Tgq 
compared to r,^ = 278(^^)1^4^ for flie Earth. 



= 270(ij^)i/% to be 



Recent 3-D atmospheric models dHeng et alJl201 la llF) sug- 
gest that G1667Cc most likely has close to 1, even in the case 
it is tidally locked and always shows the same hemisphere to 
its star, since they find that a modestly dense Earth-like atmo- 
sphere ensures a full re-distribution of the incoming energy. That 
conclusion, a fortiori, also appl ies for the presurnably denser 
atmosphere of a 4.25Me planet dWordsworth et alJlIoTll) . Wifli 



Teci ~ 270(1 - A)1/4k, G1667CC is therefore in the HZ for any 
albedo in the [0-0.83] range. 

Detailed 3-d i mensional atmospheric simulations (e.g. 
IWordsworfli et all 120111; iHeng & Vosl 1201 Ih will need to be 
tuned to the characteristics of Gl 667Cc to ascertain its possible 
climate. In the mean time, this planet, which receives from its 
star 89% of the bolometric solar flux at Earth, is a very strong 
habitable planet candidate. 



6.2. Particularity of ttie habitable zone around M dwarfs 

Two main differences between planets in HZ of solar like stars 
and M dwarfs are often pointed out. Firstly, a planet in HZ of 
an M dwarf is closer to its star, and therefore subject to more 
intense tidal forces. As a result, it is likely to quickly be cap- 
tured into a spin-orbit resonance. The Solar System demon- 
strates, however, that this does not necessarily imply that it will 
be forced into synchronous rotation. The final equilibrium ro- 
tation of a tidally influenced planet depends on both its orbital 
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Table 5. Six commensurable solutions (with similar reduced ^[)^) for the 2-keplerian orbit + 1 drift model 
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7.2002±0.0010 


4.01+0.26 


54443.33+0.53 
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c 


28.142±0.031 
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c 


90.24±0.13 
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c 
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c 
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m/s/yr 
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per degree of freedom 



Mercury, for instance, has been captured into the 3:2 , rather than 
1:1, spin-orbit resonance (* Correia & Laskaiil2004l) . and Venus 
has altogether escaped capture into a resonanc e because thermal 
atmospheric tides counteract its interior tides dCorreia & Laskaj 
I2QQ 3). Whatever the final spin-orbit ratio, the tidal forces will 
influence the night and day succession, and therefore the ch- 
mate. As discussed above however, energy redistribution by an 
atmosp here at least as dense as that of th e Earth is efficient 
JWordswordi et al.ll201ll:lHeng et al.ll2011bllah . and will prevent 
glaciation and atmospheric collapse on the night side. 

The second major difference is stellar magnetic activity, 
through its dependence on stellar mass. M dwarfs, on average, 
are much more active than solar-like star. This results from the 
compounding of two effects: lower mass stars have much longer 
braki ng times for stellar r otation dPelfosse et al J [19981 iBarnesI 
2003l: lDelorme et al.l201 Ih. and for the s ame rotation period they 
are more active (Kiraga & Stepien'2007^. As a result, a planet in 
HZ of an early M dwarf receives intense X and UV radiation for 
10 times longer than the ~ 1 00 Myrs which the solar system spen t 
close to a very active Sun (iRibas et alJl2005HSelsis et alJl2007h . 

10 



tu 0.4 





5. 10.0 20. 50. 100.0 200. 500. 1000. 



Period [days] 



Fig. 8. Periodograms for the FWHM of the correlation peak 
for the Gl 667C HARPS measurements. A clear peak at ~105 
days is well visible and is interpreted as the rotational period of 
Gl 667C. The Une show the 10% and 1% FAP level. 



eccentricity and the density of its atmosphe re ("Doyle etal.'1993; 
ICorreia et al.ll2008l: ICorreia & Laskaill20r0l: iHeller et al.ii201 L) . 
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Stellar X and UV radiation, as well as coronal mass ejections 
(CMEs), could potentially cause a major planetary atmosphere 
escape (see Scalo et al. 2007; Lammer et al. 2009, for reviews). 

Several theoretical studies have looked into the evolution 
of planetary atmospheres in M-dwarf HZ: do they lose essen- 
tial chemical species, like H2O, and can they even be com- 
pletely eroded out? The problem is complex, with many uncer- 
tain parameters, and the span of possible the answers is very 
wide. Amongst other uncertain factors, the escape ratio sensi- 
tively depends on: (1) the intensity and frequency of CMEs for 
M dwarfs (which is poorly known, and may be rnuch smaller 
that i nitialy claimed ( van den Oord & Do vld [19971 IWood et al] 
l2005h ). (2) the magnetic moment of the planet (a strong magne- 
tosp here will protect the planet against CME-induced Ion Pick 
Up dLammer et al.l l2007h ). and (3) details of the atmospheric 
chemistry and comp osition (throu gh atmosphere-flare photo- 
chemical interaction (Segura et al. 2010) and IR radiative cool- 
ing from vibrational-rotational band (Lammer et al., ,2007)). In 
two U miting cases of planets in M-dwarfs HZ. iLammer et al] 
(I2007h conclude that the atmosphere of an unmagnetized Earth- 
mass planet can be completly eroded during the Gyr-long active 
phase of its host star, while TianI (120091) find that the atmosphere 
of a 7 Earth-mass super-Earth is stable even around very active 
M dwarfs, especially if that atmosphere mostly contains CO2. 
Continuous outgassing of essential atmospheric species from 
volcanic activity can, of course, also protect an atmosphere by 
compensating its escape. 

These differences imply that a planet in the habitable zone 
of an M dwarf is unlikely to be a twin of th e Earth. Habi t abihty 
however is not restricted to Earth twins, and Ba rnes et al 1 (120101) 
conclude that "no known phenomenon completely precludes the 
habitability of terrestrial planets orbiting cool stars." A massive 
telluric planet, like G1667Cc(M2. sin/ - 4.25 Me), most likely 
has a massive planetary core, and as a consequence a stronger 
dynamo and a more active volcanism. Both factors help pro- 
tect against atmospheric escape, and super-Earths may perhaps 
be better candidates for habitability around M dwarfs than true 
Earth-mass planets. 

6.3. Gl 667Cc compared with the another known planets In 
the habitable zone 

HARPS has previously discovered two planets inside HZ. 
Gl 581d (Udry et al. 2007; Mayor et al. 2009) (M2. sin / = 7 M®) 
receives from its M3V host star just ~ 25% of the energy that the 
Earth receives from the Sun, and is thus located in t he outer hab- 
itab le zone of its star Recent detailed models ( Word sworth et al.l 
I2OII ) confirm that Gl 581d can have surface liqui d water for 
a wid e range of plausible atmospheres. HD85512b jPepe et aTl 
120 1 Ih is a M2. sin i = 3.5 planet in the inner habitable zone 
of a K5-dwarf, and receives -twice as much stellar energy as 
the Earth. It can harbor surface liquid water if it is covered by at 
lea st 50% of h ighly reflective clouds (iKaltenegger et al.ll201 ih . 

IVogt et aL (2010) announced another super-Earth in the HZ 
of Gl 581, which they found in an analysis combining HARPS 
and HIRES radial velocity data. The statistical significance 
of that detection was, however, immediately questioned (e.g. 
lTuomill20rih . and an extended HARPS dataset now demon- 
strates that the planet is unlikely to exist with the proposed pa- 
rar neters ( Forveille et al]| 201 ih . 

iBorucki et al.l (1201 l a) announced 6 planetary candidates in 
the HZ of Kepler targets with a radius below twice that of 
Earth, which they adop t as a nominal limit betw e en tell uric and 
Neptune-like planets. iKaltenegger & Sasselovl (|20I ih reduce 



this number to 3 HZ planetary candidates, by discarding 3 plan- 
ets which they find are too hot to host water if their coverage by 
reflective clouds is under 50%. Borucki et al. (201 1bi) recently 
confirmed Kepler 22b (listed as KOI-87.01 in iBorucki et al.l 
as a planet, making it the firs t planet with a measured 
radius orbiting in a habitable zone. iKaltenegger & Sasselovl 
(201 1) however discarded Kepler 22b from their candidate list 
because its radius is above 2 Earth-radii. 

To summarize, two planets with measured minimum masses 
in the range for telluric planets (Gl 58 Id and HD85512b) are 
know to orbit in habitable zones. Another one, with a measured 
radius slightly above the nominal limit for a rocky planet, orbits 
in a similar location (Kepler 22b). Finally, 3 KEPLER candidates 
with radii corresponding to telluric planets and positions in hab- 
itable zones currently await confirmation. Gl 667Cc, announced 
in Bonfils et al. (2011a) and discussed in detail here, is actually 
the most promising of those for holding conditions compatible 
with surface liquid water. Receiving ~10% less stellar energy 
than the Earth and with a minimum mass of M2. sin ; = 4.25 Mg, 
it is very likely to be a rocky planet in the middle of the habitable 
zone of its star. 

7. Summary and conclusions 

In this paper we anal ysed in detail 3 super-Earths announced in 
iBonfils et al.1 (l201Iah : Gl 433b, Gl 667Cb and Gl 667Cc. One, 
G1667Cc, is a M2.sin; - 4.25 planet in the middle of the 
habitable zone of a M1.5V star It is to date the known extra- 
solar planet with characteristics closest to Earth, but does not 
approach being an Earth twin. The main differences from Earth 
are a significantly higher mass and a different stellar environ- 
ment, which potentially can have caused divergent evolutions. 

Host stars of giant planets are preferentially metal-rich (e.g. 
Santos et al. 2001. 2004 ; Fischer & Valenti 2005) . The detection 
rate of planets with masses under 30-40M{ b, by contrast, doe s 
not clearly correlate with stellar metallicity (iMavor et al.ll201 ll) . 
This context makes the discovery of two super-Earths around 
Gl 667C ([Fe/H]~-0.6, Sect. 3) less surprising, but it remains 
one of the most metal-poor planetary host known to this date, 
with just 6 planetary host stars listed with a lower metallicity i n 
the Extrasolar Planets Encyclopaedic^ flSchneider et al.ll201 lb . 
The metallicity of Gl 433, [Fe/H]~ -0.2, is close to the median 
value of solar neighbourhood, and therefore unremarkable. 

Gl 667Cb and c orbit the outer component of a hierarchical 
triple system. Fewer than 10 other planetary systems, listed in 
iDesidera et al.l (l201 iV share that characteristics. These authors 
conclude that planets occur with similar frequencies around the 
isolated component of a triple system and around single stars. 
The planets around Gl 667C, however, are unusual in orbiting 
around the lowest-mass component of the system. 

The three planets discussed here enter in our Bonfils et al.l 
(1201 lal) statistical study, which establishes that super-Earths are 
very common around M dwarfs. Our HARPS survey of ~100 
stars has, in particular, found two super-Earths in habitables 
zones, Gl 58 Id and Gl 667Cc, even though both planets are lo- 
cated in parts of the mass-period diagram where its detection 
completeness is under 10 percent. This clearly indicates that 
supe r-Earths are common in the habitable zones of M dwarfs, 
with iBonfils et al.1 (l2011ah a 42;^jj% frequency. Future instru- 
ments optimized for planet searches around M dwarfs, like the 
SPIRou (on CFHT) and CARMENES (at Calar Alto observa- 
tory) near-IR spectrographs, will vastly increase our inventory 

' http:/; http://exoplanet .eu/index.php I 
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of such planets around nearby M dwarfs. Such instrument should 
be able to identify ~50-100 planets in the habitable zones of M- 
dwarfs, and with a 2-3% transit probability for those, to find at 
least one transiting habitable planet around a bright M dwarf. 
Such radial-v elocity educated app roach is akeady undertaken 
with HARPS jBonfils et al.l2011bl) . In this context we sti-ess that 
a photometric search for transits must be carried out shortly for 
Gl 667Cc, which have a 2.5% probability to occur. 
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Table A.l. Radial-velocity measurements and errors bars for 
Gl 433. All values are relative to the solar system barycenter, 
the secular acceleration is not substracted (see Sect. 2). 
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Table A.2. Radial-velocity measurements and errors bars for 
Gl 667C (first part). All values are relative to the solar sys- 
tem barycenter, the secular acceleration is not substracted (see 
Sect. 2). 
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Table A.3. Radial-velocity measurements and errors bars for 
Gl 667C (first part). All values are relative to the solar sys- 
tem barycenter, the secular acceleration is not substracted (see 
Sect. 2). 



Table A.4. Radial-velocity measurements and errors bars 
Gl 667C (third part). All values are relative to the solar 
tem barycenter, the secular acceleration is not substracted 
Sect. 2). 
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